Neutron-Rich Nuclei in Heaven and Earth 
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An accurately calibrated relativistic parametrization is introduced to compute the ground state 
properties of finite nuclei, their linear response, and the structure of neutron stars. While similar 
in spirit to the successful NL3 parameter set, it produces an equation of state that is considerably 
softer — both for symmetric nuclear matter and for the symmetry energy. This softening appears to 
be required for an accurate description of several collective modes having different neutron-to-proton 
ratios. Among the predictions of this model are a symmetric nuclear-matter incompressibility of 
K — 230 MeV and a neutron skin thickness in 208 Pb of R„ — R p — 0.21 fm. Further, the impact 
of such a softening on the properties of neutron stars is as follows: the model predicts a limiting 
neutron star mass of M max = 1.72M©, a radius of R = 12.66 km for a "canonical" M = 1.4M© 
neutron star, and no (nucleon) direct Urea cooling in neutrons stars with masses below M=1.3M0. 

PACS numbers: 21.10.-k,21.60.Jz,26.60.+c 



The quest for the equation of state (EOS) of neutron- 
rich matter — which is likely to lead to the discovery of 
exotic phases of matter — is an exciting problem that 
permeates over many areas of physics. While the search 
for novel phenomena has long been at the forefront of sci- 
ence, learning about "neutron-rich nuclei in heaven and 
earth" has experienced a recent revitalization due to re- 
markable advances in both terrestrial experiments and 
space observations. These developments, coupled to the 
promise of new facilities for the future of science, guaran- 
tee continuing discoveries for many years to come. Fig- 
uring prominently among the facilities for the future is 
the Rare Isotope Accelerator (RIA), a facility that by 
defining the limits of nuclear existence, will constrain 
the EOS at large neutron-proton asymmetries. In ad- 
dition, new telescopes operating at a variety of wave- 
lengths have turned neutron stars from theoretical cu- 
riosities into powerful diagnostic tools. For some recent 
excellent reviews on the relevance of the EOS on a va- 
riety of phenomena, such as the dynamics of heavy-ion 
collisions, the structure of neutron stars, and the simula- 
tion of core-collapse supernova, see Refs. 0, 0> S El an d 
references contained therein. 

Our aim in this letter is to construct an accurately 
calibrated parameter set that, while constrained only by 
the ground state properties and the linear response of 
a variety of nuclei, may still be used to predict some 
neutron-star observables. Such a successful paradigm is 
the relativistic NL3 parameter set of Lalazissis, Konig, 
and Ring |5j- The NL3 parametrization has been used 
with enormous success in the description of a variety 
of ground-state properties of spherical, deformed, and 
exotic nuclei. For some special cases, it has also been 
used successfully to compute the linear response of the 
mean-field ground state. In the particular case of the 
giant monopole resonance (GMR) in 208 Pb — the so- 
called breathing mode — the predicted distribution of 
strength is in close agreement with the experimental 
data [6j. Thus, it has come as a surprise that to re- 
produce the GMR in 208 Pb, accurately fit nonrelativis- 



tic and relativistic models predict compressional moduli 
in symmetric nuclear matter (K) that differ by about 
25%. Indeed, while nonrelativistic models predict K ~ 
220-235 MeV 0HI3, relativistic models argue for a 
significantly larger value K ~ 250-270 MeV [Bl TTol TTlj . 

In an earlier work, the density dependence of the sym- 
metry energy, which at present is poorly known, has been 
proposed as the culprit for the above discrepancy |12| . 
Since first proposed, other groups have tested this as- 
sertion reaching similar conclusions in in mm. in 

particular, in Refs. 0, it has been argued that a 
good description of the breathing mode in 208 Pb may be 
obtained using a large value of K — provided one com- 
pensates with an appropriately stiff symmetry energy. 
Thus, it was suggested that 90 Zr, a nucleus with both 
a well-developed breathing mode and a small neutron- 
proton asymmetry, should be used (rather than 208 Pb) 
to constrain K. With the advent of unprecedented exper- 
imental accuracy in the determination of the breathing 
mode in 90 Zr [6j, it now appears that the NL3 interaction 
overestimates the value of K 

El- 

Moreover, the alluded 
correlation between K and the density dependence of the 
symmetry energy serves as a telltale of a further problem 
with most relativistic parameterizations: an underesti- 
mation of the frequency of oscillations of neutrons against 
protons — the so-called isovector giant dipole resonance 
(IVGDR) — in 208 Pb El 03 • 

In this Letter we introduce a new, accurately cali- 
brated relativistic parametrization that simultaneously 
describes the GMR in 90 Zr and 208 Pb, and the IVGDR in 
208 Pb, without compromising the success in reproducing 
ground-state observables. To do so, however, two addi- 
tional coupling constants must be introduced. Without 
these additional coupling constants one can not describe 
the various modes without seriously compromising the 
quality of the fit El El- 

The effective field theoret- 
ical model is based on an interacting Lagrangian that 
provides an accurate description of finite nuclei and a 
Lorentz covariant extrapolation for the equation of state 
of dense matter. It has the following form [Til ITsl ITsj: 
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The Lagrangian density includes Yukawa couplings of 
the nucleon field to various meson fields. It includes an 
isoscalar-scalar <f> meson field and three vector fields: an 
isoscalar V M , an isovector b M , and the photon . In ad- 
dition to the Yukawa couplings, the Lagrangian is supple- 
mented by four nonlinear meson interactions. The inclu- 
sion of isoscalar meson self-interactions (via k, A, and, 
£) are used to soften the equation of state of symmet- 
ric nuclear matter, while the mixed isoscalar-isovector 
coupling (A v ) modifies the density dependence of the 
symmetry energy. While power counting suggests that 
other local meson terms may be equally important [l7| . 
their phenomenological impact has been documented to 
be small Il8l Il9j , so they will be not be considered 
any further in this study. 

A. Ground-state Properties 

Although in Refs [H EU it has been proven that the 
addition of the isoscalar-isovector coupling (A v ) is impor- 
tant for the softening of the symmetry energy, no attempt 
was made to optimize the various parameter sets. Fol- 
lowing standard practices [jlll4|. we use binding energies 
and charge radii for a variety of magic nuclei to produce 
an accurately calibrated set. We dubbed this set "FSUG- 
old" . Details about the calibration procedure will be pre- 
sented in a forthcoming publication. For now, we present 
in Table [I] a comparison between the very successful NL3 
parametrization ||, FSUGold, and, (when available) ex- 
perimental data. While the agreement with experiment 
(at the 1% level or better) is satisfactory — and this 
agreement extends all over the periodic table [2(| — a 
question immediately arises: given the success of NL3, 
why the need for another effective interaction having two 
additional parameters? The answer to this question is 
provided in the next section. 

B. Nuclear Collective Modes 

As alluded earlier, and argued in Refs. the suc- 

cess of the NL3 set in reproducing the breathing mode 
in 208 Pb is accidental; it results from a combination of 
both a stiff equation of state for symmetric nuclear mat- 
ter and a stiff symmetry energy. If true, this implies that 
NL3 should overestimate the location of the breathing 
mode in 90 Zr — a nucleus with a well-developed GMR 
strength but rather insensitive to the symmetry energy. 
Similarly, the energy of the IVGDR in 208 Pb, an observ- 



Nucleus 


Observable 


Experiment 


NL3 


FSUGold 


Ca 


B/A (MeV) 


8.55 


8.54 


8.54 




R ch (fan J 


3.45 


3.46 


3.42 




n. n — K p (rmj 




— U.UO 


— U.UO 


°Ca 


B/A (MeV) 


8.67 


8.64 


8.58 




H c h (tm) 


O.40 


Q Ad 
O.40 


Q A K 
O.40 




R n -R p (fm) 




0.23 


0.20 


yu Zr 


B/A (MeV) 


8.71 


8.69 


8.68 




Rch (fm) 


4.26 


4.26 


4.25 




R n -R P (fm) 




0.11 


0.09 


iib Sn 


B/A (MeV) 


8.52 


8.48 


8.50 




Rch (fm) 


4.63 


4.60 


4.60 




R n -Rp (fm) 




0.17 


0.13 


m Sn 


B/A (MeV) 


8.36 


8.37 


8.34 




Rch (fm) 




4.70 


4.71 




Rn-Rp (fm) 




0.35 


0.27 




B/A (MeV) 


7.87 


7.88 


7.89 




Rch (fm) 


5.50 


5.51 


5.52 




Rn-Rp (fm) 




0.28 


0.21 



TABLE I: Experimental data for the binding energy per nu- 
cleon and the charge radii for the magic nuclei used in the 
least square fitting procedure. In addition, predictions are 
displayed for the neutron skin of these nuclei. 

able sensitive to the density dependence of the symme- 
try energy, should be underestimated by NL3 (note that 
a stiff symmetry energy predicts a small symmetry en- 
ergy at the low densities of relevance to the IVGDR). In 
Table ITT1 rclativistic random phase approximation (RPA) 
results for the GMR (centroids) in 208 Pb and 90 Zr, and 
the IVGDR (peak energy) in 208 Pb are reported. These 
small-amplitude modes represent the linear response of 
the mean field ground state to a variety of probes 0, . 
Note that the FSUGold(NL3) parameter set predicts a 
compression modulus for symmetric nuclear matter of 
K = 230(271) MeV and a neutron skin in 208 Pb of 
Rn — Rp = 0.21(0.28) fm. The good agreement between 
FSUGold and experiment is due to the addition of the 
two extra parameters (£ to reduce the value of K and 
A v to soften the symmetry energy). However, it seems 
that a further softening of the symmetry energy could 
further improve the agreement with experiment. With 
the present parametrization this could not be achieved 
without compromising the quality of the fit. Thus, our 
prediction of R n — R p = 0.21 fm could be regarded as 
an upper bound. We note that the Parity Radius Ex- 
periment (PREX) at the Jefferson Laboratory is sched- 
uled to measure the neutron radius of 208 Pb accurately 
(to within 0.05 fm) and model independently via parity- 
violating electron scattering |22l |23| . This experiment 
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Nucleus 


Observable 


Experiment 


NL3 


FSUGold 


208p b 

90 Zr 

208 p b 


GMR (MeV) 
GMR (MeV) 
IVGDR (MeV) 


14.17 ±0.28 
17.89 ±0.20 
13.30 ±0.10 


14.32 
18.62 
12.70 


14.04 
17.98 
13.07 



TABLE II: Centroid energies for the breathing mode in 208 Pb 
and 90 Zr, and the peak energy for the IVGDR in 208 Pb. Ex- 
perimental data are extracted from Refs. and [2l|. 



should provide a unique observational constraint on the 
density dependence of the symmetry energy. 



C. Neutron Star Observables 

Having constructed a new accurately calibrated pa- 
rameter set, we now examine predictions for a few 
neutron-star properties. The structure of spherical neu- 
tron stars in hydrostatic equilibrium is solely determined 
by the EOS of neutron-rich matter in beta equilibrium. 
For the uniform liquid phase we assume an EOS for mat- 
ter in beta equilibrium that is composed of neutrons, pro- 
tons, electrons, and muons. Further, we assume that this 
description remains valid in the high-density interior of 
the star. Thus, transitions to exotic phases are not con- 
sidered here. 

However, at the lower densities of the inner crust the 
uniform system becomes unstable against density fluctu- 
ations. In this nonuniform region the system may con- 
sists of a variety of complex structures, collectively known 
as nuclear pasta While microscopic calculations 

of the nuclear pasta are becoming available [2(| |27], |28[ , 
it is premature to incorporate them in our calculation. 
Hence, following the procedure adopted in Ref. [2jJ, a 
simple polytropic equation of state is used to interpolate 
from the outer crust |30| to the uniform liquid. 

Results for the transition density from uniform to 
nonuniform neutron rich matter are displayed in Ta- 
ble IIIII These results are consistent with the inverse 
correlation between the neutron-skin and the transition 
density found in Ref 0] . This correlation suggests that 
models with a stiff equation of state (like NL3) predict 
a low transition density, as it is energetically unfavor- 
able to separate nuclear matter into regions of high and 
low densities. We now present results for a few neutron- 
star observables that depend critically on the equation of 
state [UEOI, namely, masses, radii, and composition. Ta- 
ble lllll includes predictions for the radius of a "canonical" 
1.4 solar-mass neutron star alongside the maximum mass 
that the EOS can support against gravitational collapse; 
beyond this value the star collapses into a black hole. 
These results were obtained by numerically integrating 
the Tolman-Oppcnheimcr-Volkoff equations. The con- 
siderable smaller radius predicted by the FSUGold model 
originates in its softer symmetry energy. The stiffer sym- 
metry energy of the NL3 set does not tolerate large cen- 
tral densities and produces stars with large radii. Note 



Observable 




r b U IjOlQ 


Pc (fm- d ) 


0.052 


0.076 


R (km) 


15.05 


12.66 


M max (M ) 


2.78 


1.72 


Pu™ (fm- 3 ) 


0.21 


0.47 


Murca(M ) 


0.84 


1.30 


AM Urca 


0.38 


0.06 



TABLE III: Predictions for a few neutron-star observables. 
The various quantities are as follows: p c is the transition den- 
sity from nonuniform to uniform neutron-rich matter matter, 
R is the radius of a 1.4 solar-mass neutron star, M max is the 
limiting mass, p Urca is the threshold density for the direct 
Urea process, Murca is the minimum mass neutron star that 
may cool down by the direct Urea process, and AMurca is the 
mass fraction of a 1.4 solar-mass neutron star that supports 
enhanced cooling by the direct Urea process. 

that the same physics that pushes neutrons out against 
surface tension in the nucleus of 208 Pb is also responsi- 
ble for pushing neutrons out in a neutron star |l8l Il9| . 
Further, while the sizable reduction in the limiting mass 
of FSUGold relative to NL3 is also due to the softening 
of the EOS, it is the softening of the EOS of symmet- 
ric nuclear matter — rather than the softening of the 
symmetry energy — that is responsible for this effect. 

We conclude with a comment on the enhanced cool- 
ing of neutrons stars. Recent observations by the Chan- 
dra and XMM-Newton observatories suggest that some 
neutron stars may cool rapidly, suggesting perhaps the 
need for some exotic component, such as condensates or 
color superconductors. Here we explore a more conser- 
vative alternative, namely, that of enhanced cooling of 
neutron stars by means of neutrino emission from nu- 
cleons in a mechanism known as the direct Urea pro- 
cess [Hill 00: 

n — > p ± e~ + D e , (2a) 
e~ ± p — > n ± v e . (2b) 

This mechanism is not exotic as it only relies on protons, 
neutrons, electrons, and muons — standard constituents 
of dense matter. However, it requires a large proton frac- 
tion Y p for the momentum to be conserved in the above 
reactions. As a large proton fraction requires a stiff sym- 
metry energy, it is interesting to determine if the newly 
proposed EOS is able to support such a large proton frac- 
tion. Note that in order for the direct Urea process to op- 
erate, the proton fraction must exceed 1^ = 0.111 for the 
low-density (muonless) case, and 1^ = 0.148 for the high- 
density case (with equal number of electrons and muons) . 
In Table lull we list the threshold density (p Urca ) and min- 
imum mass (Mrjrca) required for the onset of the direct 
Urea process. We note that in spite of its softer symme- 
try energy, FSUGold predicts that a 1.4 solar- mass neu- 
tron star may cool down by the direct Urea process. For 
completeness, the mass fraction that supports enhanced 
cooling in such a neutron star is listed as AMijrca- 
In conclusion: 
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1. A new accurately calibrated relativistic model 
( "FSUGold" ) has been fitted to the binding energies and 
charge radii of a variety of magic nuclei. The model is 
as successful as the NL3 set — used here as a success- 
ful paradigm — in reproducing the ground-state prop- 
erties of a variety of nuclei. Symmetric nuclear matter 
saturates at a Fermi momentum of ftp = 1-30 fm -1 (corre- 
sponding to a baryon density of 0.15 fm -3 ) with a binding 
energy per nucleon of B/A — —16.30 MeV. Yet relative 
to NL3, FSUGold contains two additional parameters, 
whose main virtue is the softening of both the EOS of 
symmetric matter and the symmetry energy. 

2. While the need for the two additional parameters 
is not manifest in the description of ground-state prop- 
erties, it becomes essential for reproducing a few nuclear 
collective modes. Specifically, the breathing mode in 90 Zr 
is sensitive to the softening of symmetric matter, the 
isovector giant dipole resonance in 208 Pb to the softening 
of the symmetry energy, and the breathing mode in 208 Pb 
to both. Incorporating these additional constraints yields 
a nuclear-matter incompressibility of K — 230 MeV and 
a neutron skin thickness in 208 Pb of R n — i? p = 0.21 fm. 

3. While the description of the various collective modes 
imposes additional constraints on the EOS at densities 
slightly above and below saturation density, the high- 



density component of the EOS remains largely uncon- 
strained. We made no attempts to constrain the EOS 
at the supranuclear densities of relevance to neutron-star 
physics. Rather, we simply explored the consequences 
of the new parametrization on a variety of neutron star 
observables, namely, masses, radii, and composition. In 
particular, we found a limiting mass of M max = 1.72M©, 
a radius of R= 12.66 km for a canonical M = 1.4M© neu- 
tron star, and no direct Urea cooling in neutrons stars 
with masses below M =1.3Mq. While the consequences 
of these results will be fully explored in a forthcoming 
publication, it is interesting to note that recent observa- 
tions of pulsar-white dwarf binaries at the Arecibo ob- 
servatory may place important constraints on our EOS. 
Indeed, Nice, Splaver, and Stairs have inferred a pulsar 
mass for PSRJ0751+1807 of M = 2.1±%£M e at a 95% 
confidence level [3^ . If this observation could be refined, 
not only would it rule out the high-density behavior of 
this (and many other) EOS, but it could provide us with 
a precious boost in our quest for the equation of state. 

We acknowledge many useful discussions with C.J. 
Horowitz and D. Page. This work was supported in part 
by DOE grant DE-FG05-92ER40750. 
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